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Abstract 

In a first paper we proposed a continuum model for the dynamics of systems of self 
propelling particles with kinematic constraints on the velocities. The model aims to 
be analogous to a discrete algorithm used in works by T. Vicsek et al. In this paper 
we prove that the only types of the stationary planar solutions in the model are 
either of translational or axial symmetry of the flow. Within the proposed model 
we differentiate between finite and infinite flocking behavior by the finiteness of the 
kinetic energy functional. 
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1 Introduction 



The dynamics of systems of particles subjected to nonpotential interactions 
remains poorly understood. The absence of a Hamiltonian for such systems, 
which generally are far from equilibrium, hampers applying the machinery of 
statistical mechanics based on the Liouville equation. Many attempts have 
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been made to investigate these systems using discrete algorithms to model 
this behavior. In nature there are many examples of such systems [1]. Since 
the discrete algorithms are hard to describe analytically it is natural also to 
consider continuum models of a hydrodynamic type. In standard hydrodynam- 
ics the relation between microscopic kinetics (Boltzmann-type equations) and 
Navicr-Stokes equations is a standard topic of research [2]. For the systems 
of interest the construction of corresponding kinetic equations based on the 
specific dynamic rules and their connection with the hydrodynamics equations 
seems to be unknown so far and is worth studying. One can expect that the 
continuous description of the collective behavior like swarming and flocking 
leads to quite unusual hydrodynamics. 

In our first paper [3] we proposed a hydrodynamic model which can be con- 
sidered to be the continuum analogue of the discrete dynamic automaton 
proposed by Vicsek et al. [4] for a system of self propelling particles. It uses 
the continuity equation 

^^^ + div(n(r,t)v(r,t)) = 0, (1) 
which implies that the total number of particles 

J n{r,t)dr (2) 
is constant. The kinetic energy of a co-moving volume 

T^^ln{r,ty{r,t)dr (3) 

is also conserved 

Using Eq. (1) and Eq. (4) it can be shown that a field u) exists such that the 
Eulerian velocity, v(r,t), satisfies: 

^v(r,t)=a;(r,t)xv(r,t). (5) 

This equation can be considered as the continuous analogue of the conservative 
dynamic rule used by Vicsek et al. [4]. 

We proposed the following "minimal" model for the field of the angular velocity 
oj (r, t) which is linear in spatial gradients of the fields n(r, t) or v(r, t): 

u (r, t) = j Ki{r- r') n(r', t) rot v(r', t) dr'+ J K2{r- r') Vn (r', t) x v(r', t)dr' . 

(6) 

uj has the proper pseudovector character. The averaging kernels Ki (r — r') 
and K2 (r — r') should naturally decrease with the distance in reahstic models. 
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They sample the density and the velocity around r in order to determine 
uj{r,t). In the first paper we concentrated on Ki. The detailed derivation 
of the above equations from the discrete models based on the automaton 
proposed by Vicsek et al. [4,5] will be the subject of a future paper. 

Note that the models based on Eqs. (l)-(6) allow solutions of uniform motion 
in the form of a solitary packet: 

n(r,i) = no(r - vot) (7) 

with vo independent of position and time. The contribution to uj due to Ki 
is zero for an arbitrary density distribution uq. The contribution due to K2 
is zero for density distributions no which only depend on the position in the 
Vo direction. In this second case it follows from the continuity equation that 
no should be everywhere constant. The density distribution Uq should be cho- 
sen such that the number of particles, and, correspondingly, the total kinetic 
energy are finite. The solutions of such type also were found analytically in 
[6] and observed in simulations [7] . Note that such solutions exist not only in 
nonlocal case like in [6] but also for the local model which we consider below. 

Within the first order of perturbation theory on small deviation of density and 
velocity fields the sofitary solution given by Eq. (7) shows neutral stability; 
i.e. the perturbations grows linearly for small t. 

We restrict our discussion to the simple case of averaging kernels, which are 
5-functions: 

Kj (r - r') = Sj S{r- r') , where j = 1 or 2. (8) 

We will call this the local hydrodynamic model (LHM). In the first paper, 
where we only considered Ki, we scaled Ki by dividing by and the den- 
sity n by multiplying with |si|. This made it then possible to restrict the 
discussion to Si is plus or minus one. The disadvantage of this scaling pro- 
cedure is that it changes the dimensionality of Kj and n. For two kernels it 
becomes impractical. We note that Sj is given by: 

Sj = / Kj (r) dr. (9) 

For the local model Eq.(6) reduces to: 

(jj(r,t) = Sin(r,t)rotv(r,t) + S2Vn(r,t) x v(r,t). (10) 
and Eq. (5) for the velocity becomes 

(r, t) = si n(r, t) rot v(r, t) x v (r, t) + S2 (Vn (r, t) x v(r, t)) x v (r, t) . 

(11) 
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Note that the second term on the right-hand- side of Eq. (11) corresponds to 
the rotor chemotaxis force when the number density is proportional to the 
attractant density introduced in [5]. 

In the following section we will show that the only stationary solutions for 
LHM with ijj field given by Eq. (10) are either the solutions of uniform motion 
(see Eq. (7)) or the radially symmetric planar solution which will be considered 
in detail in the following section. 

In the second section we investigate the properties of the stationary radially 
symmetric solutions of the local hydrodynamic model for some special cases. 
Conclusions are given in the last section. 



2 Possible types of stationary states for the local hydro dynamical 
model 

The equations of motion to be solved are Eqs. (1) and (11). In order to find 
a class of 2D stationary solutions we consider this problem in a generalized 
curvilinear orthogonal coordinate system (m,v), which can be obtained from 
the Cartesian one (x, y) by some conformal transformation of the following 
form: 

u^iv = F{z), (12) 

where F(^z) is an arbitrary analytical function oi z = x ^ iy. In a curvilinear 
orthogonal coordinate system the fundamental tensor has a diagonal form, 
9ik = 9ii^ik^ where the indices i,j are either u ox v. The square of linear 
element in conformal coordinates is 

ds'^ / A du^^dv'^, (13) 
D{u,v) ^ ' 

where 

^, > diu.v) ( du\'^ f du\^ f dv\^ ( dv\^ 1 ,^ 

^("■") -wd-r[di)^ i%j - [si] ^[ii)-T9 '''' 

is the Jacobian of the inverse conformal transformation from the arbitrary 
curvilinear orthogonal to Cartesian coordinates, {u,v) — > {x,y). Furthermore 
g = Quugw is the determinant of the metric tensor. For a conformal transfor- 
mation = Qy^u = 

The differential operations are given by the following expressions [8]: 
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\ ou ov J 



div a — D 
rota — D 



d ( a„ 



dv 



d I _a^\ d ^ an 
du[,/D dv 



D. 



d^4> ^ d\ 



(15) 
(16) 
(17) 
(18) 



Here and are orthonormal base vectors in the directions of increasing u 
and V respectively. These base vectors are functions of the coordinates u and 
V. The projections of the vcctorfield a on these directions are a^j = a ■ and 
= a - e^. Using Eqs. (15)-(18) for the velocity field given by v = v„ e„+v^ e^, 
one obtains: 



(v • V) V = 



divuBu) , 5(v„e„) d{vyey) 5(v^e^ 
v«, — h v^ — h Vu — h v„- 



du 



dv 



du 



dv 



rot V X v — D 



d f Vv \ d f Vu 



(Vn X v) X V ^Vd 
div(nv) = D 



du \VD J dv \^/D 
dn dn 



ou ov 



(-v^e„ + v„e^), 



d I nv„, \ ^ d I nvr, 



du \y/D J dv \^/D, 



0. 



(19) 

(20) 
(21) 

(22) 



Substituting Eqs. (19)-(22) into Eqs. (1) and (11) we obtain the following 
system of equations, which determines all possible stationary flows for the 
LHM: 

dvu , ^ dvu , f J. , sindlnD dn\ 2/ SiudlnD dn\ dy^ 

v.— +(l-.,n)v.— +V.V. (^/s + — ^ - s,^j+v^ (^/4 - — ^ + S2g^j+s,ny.^ 

(23) 



dvy dvy ( SindlnD dn\ 2 f r sindlnD dn\ 9v„ 

v.— +(l-.,n)v.— +v„v. [h + ^^-^^Q^) +v„ [h - + ^^9^) +^^^^« - ' 

(24) 

d ( nv,, \ d ( riY,,, \ , , 



du \ / dv\ -/d 

where 
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fi{u,v) = 



du 
du 



f2{u,v) = 

Mu,v) = 



dv 
dv 



(26) 



Now let us consider the case of " coordinate flows" , when the flow is directed 
along one of the families of coordinate lines u, v for example along w-coordinate 
lines and is given by v = (0, v^(u, v),0), the density distribution isn — n{u, v). 
The case of a velocity field v = Vu^u is equivalent (just interchange u and v). 
Prom Eq. (23) we have: 



where 



I{u,v) 



= C exp 

1 d \nD 
2 



j I{u,v)du , 

/4 828 Inn 



du Sin 
Equations (24) and (25) take the form: 



si du 



(27) 



(28) 



dvv_ 
dv 



0, 



d I ny 
dv 

and lead to 



D. 



n(u,v) — h(u)^J D(u, v), 
where h{u) is an arbitrary function of u. Taking into account that: 



h{u,v) = 
and Eq. (31) we obtain: 



dv 



Dd{l/D) _ ld\nD 
~2du' 



I{u,v) 



1 d \nD 

2 du 



Sih{u)vD si 
Note that as it follows from Eq. (29) 

For the integrand in Eq.(27) this implies that 

I{u, v) = I{u). 



du 



S2 ( ^d lnh{u) ( d \nD 



(29) 
(30) 

(31) 
(32) 



-1 



du 



du 



+ 1 



(33) 
(34) 

(35) 



Therefore from Eqs. (33) and (35) we can conclude that the function D[u,v), 
which determines the coordinate system, depends only on one variable, D ~ 
Diu). 
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In the case of conformal coordinates, defined by the metrics in Eq. (13), the 
Gaussian curvature of the surface is given by [9] : 

K = ^DAlnD. (36) 

For a planar flows the condition K — leads to the following 

AlnD = 0. (37) 

Using the expression for the Laplacian in the conformal representation (see 
Eqs. (18)) and taking into account a fact that D — D{u) one finds for (37): 

D = exp [ciu + C2] , (38) 

where ci^2 are arbitrary constants. The case Ci = determines a Cartesian 
coordinate system, which related to a linear class of stationary flow. The case 
ci 7^ determines a polar coordinate system [8], which corresponds to an 
axially symmetric (or vortical) type of flow. 

Finally the velocity field for the LHM, with si takes the form: 
d InD 



Yv (u) = C exp 



du 



1 L__faf2^i^f^^V' + i| 

Sih{u)VD si \ du \ du J 



(39) 

Thus it is proved that for the case si 7^ the only stationary solutions are 
those either with planar or axial symmetry of the fiow. 

The case Si = (the LM2) is specific because as it follows from Eqs. (23) and 
(24) the velocity field Vv{u) is arbitrary while the density takes the form: 

n^-J—]nD + no. (40) 

2 S2 

The statement about the symmetry of the stationary solutions for such a model 
is the same as that proved above for the case Si 7^ 0. Note that the parameter 
A = S2/S1 can be considered as the weight factor of the rotor chemotaxis 
contribution. 



3 The properties of radially symmetric stationary solutions 



In this section we investigate the stationary radially symmetrical solutions 
for different cases of the local hydrodynamical models. In our first paper we 
considered the S2 = case, which we called the local model one (LMl). Other 
models correspond to cases si = and si — S2 — s which we will call the 
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LM2 and the LM12 respectively. We consider the finite and infinite fiocking 
stationary states for these models. It is natural to differentiate between these 
two cases by the finiteness of two integrals of motion - the total number of 
particles Eq. (2) and the kinetic energy Eq. (4). The infinite flocking is asso- 
ciated with N infinite but finite T, while finite fiocking naturally corresponds 
to both N and T finite. Note that in the finite fiocking behavior one may 
consider two cases with respect to the compactness of the n(r). Compactness 
means that the density has some upper cut-off beyond which it can be put 
zero. 

3. 1 The properties of the stationary solutions of LMl 

In our previous paper [3] we considered the stationary solutions in LMl and 
obtained the for velocity field profile 



Here r is the radial coordinate and Lp is the angular one, Tq is some radius, 
which for vortex-like solutions plays the role of lower cut-off radius of the 
vortex and determines its core. 

For vortex-like solutions the constant Ci in Eq.(41) is determined by the 
circulation of the core 



The spatial character of the solution given in Eq.(41) strongly depends on the 
sign of the parameter si. The finiteness of integrals of motion Eq. (2) and 
Eq. (3) is guaranteed by either the fast enough decrease of the density n(r) 
at r — >^ oo or its compactness (n(r) as a function has finite support). 

Let us consider the finite fiocking behavior (FEB), which is characterized by 
both and T finite. If at r — cxd asymptotically n{r) ~ r~° , where a > 2 
the total number of particles is finite. Then at such a behavior of n(r) the 
total kinetic energy is finite only if si < 0. 

In a case Si > the condition of finiteness for the kinetic energy and the total 
number of particles is fulfilled only if n{r) has finite support. As an example 
we may give: 



where R is the upper cut-off radius. Substituting Eq.(43) into Eq.(41) one 




(41) 




(42) 




(43) 
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obtains: 




arctanh ~ 



ro 



(44) 



The corresponding profiles for the velocity v = v^^e^ at difi^erent ratio R/tq 
are shown in Fig. 1. Note that for the case considered at R/tq > 2 we get 
the monotonic profiles of the velocity which are similar to those observed in 
experiments [5]. 




Fig. 1. The velocity profiles for V* = 2nRv^{r /R) in the LMl for different ratios 
R/ro and si^i = 1. 



The infinite flocking behavior (IFB) is characterized by the N inflnite and T 
finite. For the case Si > no physical solutions exist with such a behavior. For 

the case Si < slowly decaying density distributions n(r) oc r^" at r — > cxd 
with < « < 2 are consistent with the finiteness of T. These statements are 
summed up in Table 1. 



LMl 


si > 


si < 


FFB (AT < oo, T < oo) 


compact support 


a > 2, no compact support 


IFB (iV = oo, T < oo) 


no physical solutions 


< a < 2 



Table 1 



Properties of the stationary solutions of the LMl. 



3.2 The properties of the stationary solutions of LM2 



For the case that Si = and S2 finite one may also construct a radially 
symmetric stationary planar solution. In polar coordinates from Eq. (40) we 
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get: 

1 T 

nir) = — In — , (45) 

and one can choose the velocity field V(p(r) arbitrarily. For positive values of 
S2 this density is positive for r > ro and for negative values of S2 it is positive 
for r < ro- So for positive values of S2 the density profile becomes 

n(r) = —e{R-r) In— (46) 

and for negative values of S2 it becomes 

n{r) = - — j- 9 (ro — r) In — . (47) 

1^21 



The results about finite and infinite flocking behavior are in Table 2. 



LM2 


si > 


si < 


FFB (,Y < 3c, oc) 


uo physicril solulion 


compaci sui)p()ri 


IFB (AT = oo, T < oo) 


no physical solution 


no physical solution 



Tabic 2 



Properties of the stationary solutions of the LM2. 



3.3 The properties of the stationary solutions of LM12 



The third case which is expedient to consider is Si = S2 = s- In that case, 
according to Eq. (10), the field uj is coupled to the number density fiux j = n v: 



u? (r, t) — s rot j 
so that Eq. (11) for the velocity is 

^v(r,i) = s rotj x v (r,t) . 

For a radially symmetric stationary planar solution this gives 

d 



with 



yjr) 



C 



2'Kr n{r) 



exp 



s Jro r n (r j 



(48) 
(49) 

(50) 



(51) 



as a solution. The constant C is determined by the circulation of the core 

C 



vdl 



r=ro 



n{ro) 



(52) 
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The properties of finite and infinite flocking behavior for this model are the 
same as those for the LMl (see Table 1). 



4 Conclusions 

In this paper we consider the properties of the stationary 2D solutions of 
the LHM proposed in [3]. We established that the only possible stationary 
solutions in the model are those with translational or axial symmetry. The 
cases of finite and infinite flocking behavior are considered for different specific 
types of the LHM. It is shown that the case Si = (LM2) is specific in a sense 
that there is only one density distribution, for which many velocity profiles 
can be realized. In general case (si ^ 0) one is free to choose axially symmetric 
density distribution which the velocity profile depends on (Eqs. (31) and (39)). 
Note that in this respect the general case is similar to the LMl considered 
earlier. 
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